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Abstract 

In this paper, we construct under general assumptions the stochastic dynamics of an inter¬ 
acting particle system in a bounded domain 0 with sticky boundary. Under appropriate 
conditions on the interaction the constructed process solves the underlying SDE for every 
starting point in the state space. Moreover, we also obtain a solution for q.e. starting point 
in the case of singular interactions which generalizes former results. Finally, the setting is 
applied to the case of particles diffusing in a chromatography tube. 
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1 Introduction 


In |FGV14| a sticky reflected distorted Brownian motion on [0, oo) n , n € N, is constructed via 
Dirichlet forms and applied to stochastic interface models. Afterwards, the connection to random 
time changes and Girsanov transformations is investigated in |GV14b| . In particular, strong Feller 
properties of the transition semigroup of a process associated to the underlying Dirichlet form 
are proven such that the existence result for weak solution of the underlying SDE of |GV14b| is 
improved under appropiate assumptions on the drift. Moreover, the Dirichlet form construction of 
sticky reflected distorted Brownian motion on the half-line [0, oo) is generalized to general bounded 

domains O C M d , d G N, in | GV14a] such that the setting even allows a diffusion on dQ. In the 

- N 

present paper, we construct and analyze diffusions on O , N £ N, with sojourn on the boundary. 
This type of diffusion desribes naturally a system of interacting particles with sticky boundary. 
In the independent case, i.e., the case without interaction, the setting reduces to N independent 
diffusions, where each diffusion is of the type considered in |GV14b] , We define the corresponding 
Dirichlet form and present the connections to random time changes and Girsanov transformations. 
Moreover, we calculate the corresponding L 2 -generator for smooth functions and establish in this 
way the connection to the underlying martingale problem and SDE. 

The construction allows very weak assumptions on the interaction and moreover, illustrates some 
effects which do not appear in the case of interacting particle systems with absorbing or reflecting 
boundary conditions. For example, in the case N — 1 without drift the invariant measure for the 
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2 Preliminaries 


sticky reflected Brownian motion on fl is given by A + cy where A denotes the Lebesgue measure 
on fl and a the surface measure on dPl. Then, the invariant measure for general N G N with 
additional drift is given by 

N 

p = + ^ 

2—1 

where g is a suitable density. In the case of absorbing or reflecting boundary conditions the in¬ 
variant measure can be derived similarly, but the surface measure does not appear. Hence, the 
structure of the product measure is much simpler. In this case, one obtains the Lebesgue measure 
on f2 := fl N and previous results apply. The only problem is that usually the boundary regularity 
decreases, since II possesses corners. In the case of the sticky boundary condition it is not possible 
anymore to reduce the setting for general IV G N to the case N — 1. Therefore, it is necessary to 
analyze the structure of the problem in detail. 

The investigated system of SDEs is of the form 

( ;x; =i a (xi)(dB' t + i(L^(x;) + ^(x«))a) - ir(x;)^(x;> n(x\)dt 
+s lr(X<)( dB, ’■* + (x;) + Z**C*.))*). * = 1, • • •. N 

dB ,P = P(Xj) o dB\ 

X 0 =x G if, 

where 5 G {0,1}, (B t ) t > 0 , B t = (B ],..., B^), is an A r d-dimensional standard Brownian motion, n 
is the outward normal vector and P is the projection on the tangent space. The particle interaction 
is given by V* In (j) and Vr,* In (j), i — 1,..., N, where Vp denotes the surface gradient. The precise 
definitions of V* and Vr,i as well as n and P are given in Section [?T1 The densities a, and /?j, 

1 = 1 ,N, are only assumed to be continuous and to fulfill a weak differentiability condition 
whereas (j) is C 1 (see also Condition 14.21 and Theorem 14.101) . Note that the drift is nevertheless not 
necessarily Lipschitz continuous, since the densities are allowed to vanish on a set of measure zero. 
A similar system of SDEs has been investigated in | IGra88] and applied to a model for molecules 
diffusing in a chromatography tube. We also consider such kind of applications and extend previ¬ 
ous results to the case of singular interactions. 

Our paper is organized as follows: In Section [2] basic notations are explained and some previous 
results are stated. In Section [3] the underlying Dirichlet form is constructed and afterwards, 
in Section [T] the associated diffusion is analyzed and the relations to random time changes and 
Girsanov transformations are presented. Finally, we apply the results in Section 0 

2 Preliminaries 

2.1 General notation 

Throughout this paper, fl C d > 1, denotes a nonempty bounded domain such that its bound¬ 
ary T := f)Q is of Lebesgue measure zero. In the case S = 1 we assume that d > 2. The standard 
scalar product in R n , n G N, is given by (•, •) and norms in KG by | -1 (in particular, for the modulus 
in R; eventually labeled by a lower index in order to distinguish norms). Similarly, || • || denotes 
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2 Preliminaries 


norms in function spaces. The metric on R d induced by the euclidean metric is denoted by d euc . 

For a vector x E fl N , JV e N, we use the representation x = (x 1 ,... ,x N ), where x l E 0, i = 
1,..., N, is represented in the form x 1 = (x\, ..., x l d ). We denote by V the gradient of a smooth 
function and by d x %, i = 1,... ,N, k = 1 ,,d, its partial derivatives. In the case JV = 1 we 
simply write 3*. for k — 1,..., d. By V,;, % = 1,... ,N, we denote the d-dimensional vector given 
by the partial derivatives with respect to the coordinates x\, k — 1,..., d. Moreover, V 2 denotes 
the Hessian for functions mapping from subsets of R d to R and A = Tr(V 2 ) the Laplacian. V 2 
and Aj are dehned analogously. In the case of Sobolev functions we use the same notations in the 
weak sense. 


2.2 Submanifolds in the euclidean space 

In the following, the boundary T of Q is said to be Lipschitz continuous (respecktively C' fc -smooth) 
if Definition 2.1 of |GV14a] holds, i.e., F is Lipschitz continuous (respectively C'Wsmooth) if D 
is locally below the graph of a Lipschitz continuous (respectively C k -) function and the graph 
coincides with T. In this case, we also simply say that F is Lipschitz (respectively C k ) or that 
fl has Lipschitz boundary (respectively C k -boundary). Moreover, the surface measure on F is 
denoted by a and the (outward) normal vector at a point x E F is denoted by n(x) (supposed the 
boundary is smooth at x). 

Remark 2.1. The definition of n can be extended to a neighborhood of x and n is differentiable 
near x if T is C 2 . 


Definition 2.2. Let x E F be such that the outward normal n(x) exists. Define 

P(x) := E — n(x)n(xy E R dxd , 

where E is the d x d identity matrix. We call P(x) the orthogonal projection on the tangent 
space at x. Note that P(x)z = z — (n(x), z ) n(x) for z E R d . 

Definition 2.3. Let f E C l (Q) and x E L. Then we define (whenever T is sufficiently smooth at 
x) the gradient of / at x along T by 

V r f(x) := P(x)Vf(x) 

and if f E C 2 (Q) the Laplace-Beltrami of / at x by 

A r f(x) := Tr(V 2 /(x)) = div r V r /(x) = Tr (P{x)V(P{x)Vf(x))), 

where divr$ := Tr(PV$) for = ($i,. . ., E C' 1 (h2; R d ) with = J<I> = (V$i| .. . |V$d). 
Analogously, we define higher derivatives of order k E N. In this way, let C fc (ro) be the space of 
continuously differentiable functions on To obtained by restriction of C k (Fl) -functions, where T 0 is 
an open subset of T in the subspace topology. As usual, set C' oo (T 0 ) := C^To). Moreover, 
in the case that n is differentiable at x we define the mean curvature of T at x by 


k(x) := divr n(x). 

Remark 2.4. For smooth functions, we have the divergence theorem 

J (<f>, Vr<?) da = — J div r <l> g da , 

where 4> is Revalued ( se e e.g. |Tayll Chap. 2, Proposition 2.2]). 


( 2 . 1 ) 
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The following lemma follows easily by calculation: 

Lemma 2.5. Assume that T is C 2 -smooth. Then 

(PV/.P = -Kn. 

Definition 2.6. Let To be an open subset of T in the subspace topology. The Sobolev space 

H L ’ k (To), k> 1, is defined by C' 1 (T 0 )" ^ Hl ’ fe(r<)) C L k (T 0 ;a), i.e., the closure C' 1 (T 0 ) with respect to 
the norm 

II • lltf^To) := (II • \\ k L^(To-,a) + II V r • llz,fc(r 0 ;o-)) * • 

Remark 2.7. id 1 ,fc (r 0 ) can also be charaterized as the space of functions which are in local coordi¬ 
nates in the corresponding Sobolev space. 

If / G H 1,k ( To) and (/ n )neN is an approximating sequence of smooth functions, Cauchy in H 1,k (T o), 
we call the L k (T 0 ; cr)-limit of (Vr/ n )neN the weak gradient of / and denote it by Vr/- In the case 
T 0 = T, (12.ip transfers from f n to / using a continuity argument provided that d> G L k (T; er) for 


2.3 Brownian motion on manifolds 

We shortly recall some facts about Brownian motion on T. For details about stochastic analysis 
on manifolds, we refer to |HT94] . (H su 02] and |IW89] . 

By definition, Brownian motion (P/) t > 0 on T is a T-valued stochastic process that is generated 
by |Ap, in analogy to Brownian motion on M. d , in the sense that (Bf ) t >o solves the martingale 
problem for (|Ar, C°°(T)). We recall the following: 

Lemma 2.8. Let T be (7 2 -smooth. Then a solution of the Stratonovich SDE 

dx t = p(x t ) o dB u x 0 g r, 

is a Brownian motion on T, where (B t ) t >o is a Brownian motion in M. d . 

Proof. See |Hsu021 Chap. 3, Sect. 2], □ 

Remark 2.9. Note that the dimension of the driving Brownian motion (B t ) t > 0 is strictly larger 
than the dimension of the submanifold T and hence, according to |Hsu02] the driving Brownian 
motion contains some extra information beyond what is usually provided by a Brownian motion 
on T. Furthermore, a solution of the above SDE is naturally T-valued, since P(x)z is tangential 
to T at x for every x G T and z G M/ In our application, it is natural to construct a Brownian 
motion on T by means of a d-dimensional Brownian motion, since a Brownian motion on W l is 
involved anyway. 

We also need the following result: 

Lemma 2.10 (Ito-Stratonovich transformation rule). Consider the Stratonovich integral in given 

by 

5(X t ) o dB t , 

where B = (B t ) t >o is a d-dimensional Brownian motion and S : i —> R dxd is C Yl -smooth and 

symmetric. Then the Ito form reads 

S(X t )dB t + i(( SWyS)(X t )dt . (2.2) 
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Note that in the case S = P the (12.211 can be represented in the form 

P(X t )dB t - l«n(X t )dt 

in view of Lemma 12.51 

2.4 Sticky reflected diffusions on Q 

In the following we recall the main results of [ GV14aJ . 

Assume that T = dfl is Lipschitz continuous. Moreover, assume a G L 1 (h2; A), a > 0 A-a.e., and 
(3 E L l (T] a), [3 > 0 er-a.e.. 

Define 


Q In ol + lr f3 


(2.3) 


as well as 

/i := q (A + a) = a\ + /3a. 

Note that the condition a E L 1 (D; A), a > 0 A-a.e., and /? G iAfT; a), (3 > 0 cr-a.e. is equivalent 
to q G L 1 (h2; A + a), g > 0 (A + er)-a.e.. 

Let the symmetric and positive definite bilinear form (S,T>) be given by 


S(f,g)-=T / (V/.Vs) ad\ + - (Vrf.Vrg)/3da for f,g £ V := C\Sl), 


(2.4) 


in 


where (•, •) denotes the euclidean scalar product in M. d and 5 G {0,1}. In addition, let 


as well as 


Sn(Lg) ■■= - / (v/, X7g) adX for f,gEV n := 

2 . Cl 


£r{f,g) '=2 ( Vr /’ Vr 5 f ) P da for /:.9 e v r ■= 


Note that e([D) = e(X>n) — Dr, where e : C' 1 (h2) —Y C ,1 (T) is dehned by the restriction of functions 
to T. In this terms, for f,gEV we get 

S(f,g) = S n (f,g) + 5S T (f,g). 

In order to prove closability of we need an additional assumption on the density g. Define 


R a (Pl) := {x E Pi : a l d\ < oo for some e > 0} 

J {yen-.\x-y\<e} 

and analogously Rp(T) with PI replaced by L and A replaced by a. We assume that a = 0 A-a.e. 
on Pl\R a (Pl) and additionally [3 — 0 a- a.e. on r\i? ( g(r) if 5 = l(Hamza condition). 

Under these assumptions the following holds true: 

Theorem 2.11. The symmetric and positive definite bilinear form (£,D) is denesly defined and 
closable on L 2 (0; /i). Its closure (£, D{£)) is a recurrent, strongly local, regular, symmetric Dirich- 
let form on L 2 (D; /i) . 
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2 Preliminaries 


As an immediate consequence we obtain the following theorem: 

Theorem 2.12. There exists a conservative diffusion process (i.e. a strong Markov process with 
continuous sample paths and infinite life time) 

M := (12, T , ( Tt)t>o, (X*)t>o, (0t)t>o, (P 

with state space 12 which is properly associated with (£,D(£)), i.e., for all (p-versions of) f G 
Bbifil) C L 2 (12; p) and all t > 0 the function 

09^ ptf(x) := E x (f{X t )) := / f(X t )dP x G R 

is a quasi continuous version of T t f. M is up to p-equivalence unique. In particular, M is p- 
symmetric, i.e., 

Ptf g dp = / / p t g dp for all f,gE 23&(12) and all t > 0, 

J o 

and has p as invariant measure, i.e., 



Ptf dp 


f dp for all f G 23&(12) and all t > 0. 

Jn 


If we assume the stronger conditions that T is (7 2 -smooth and a, (5 G C(12), a > 0 A-a.e. on 12, 
ft > 0 a-a.e. on F such that y/a G // 1,2 (12) and additionally yffi G H 1,2 (T) if 8 — 1, it is possible 
to determine the generator of (£, D(£)) for functions in (7 2 (12). The explicit representation of the 
generator allows to analyze the dynamics of M: 


Theorem 2.13. M is a solution to the SDE 

dX t =l n (X t )(dB t + l — (X t )dt) - lr(X f )^(X t ) n{X t )dt 
V 2 a / p 

+ 5 lr(X 1 )(dB, r + i^M(X t )*), (2,5) 

dBf — B(X t ) o dB t , 

Xq .1 , 

for quasi every starting point x G 12, where (B t ) t > o is a d-dimensional standard Brownian motion, 
i.e., 


X t = x+ [ t l n (X a )dB 8 + I 1 n (X s )l — (X s )ds 
Jo Jo 2 a 

+ 6 I l r (X.)P(X,)dB, - 5 [ l r (X»)i K (X»)n(X,)ds 


( 2 , 6 ) 


^0 

+ S J) 


>0 


[ ~(X s )l r (X s )n(X s )ds 
'o 2 p 


almost surely under P x for quasi every x G 12. 
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3 The Dirichlet form and the associated Markov process 


If we suppose addtionally to the assumptions of Theorem 12.131 that there exists p > 2 with p > | 
such that 

In I a I + $ lr ^ ^ G tf oc (Q D {,o > 0}; p) 

and cap£-({^ = 0}) = 0, we even obtain a stronger version: 

Theorem 2.14. There exists a conservative diffusion process 

M = (12, T, (fF t )t> o, (X t )t> 0 , (0t)t>o, (P®)a. e nn{ e >o}) 

with state space 12 D > 0} such that M solves \2.5\l for every xGl2fl{£>>0}. Moreover, its 
Dirichlet form is given by (£, £)(£)) on L 2 (12 D {p > 0}; p) and the transition semigroup (pt)t>o 
of M is C p -strong Feller, i.e., p t (C p (Fl fl {g > 0}; p)) C (7(12 D {p > 0}). In particular, (p t )t >o it 
strong Feller, i.e., pt(Bb(Fl fl {p > 0})) C (7(12 D {£> > 0}). Furhtermore, M has a sticky boundary 
behavior, i.e., 

lim - / l r (X s )ds > 0 

t^oo t J o 

P x -a.s. for every x G 12 fl {p > 0} such that x is in a component of Ft fl {£> > 0} intersecting T. 


3 The Dirichlet form and the associated Markov process 


3.1 General setting 

Assume that T := <912 is Lipschitz continuous. Let (G,D(Q)) be the recurrent, strongly local, 
regular, symmetric Dirichlet form on L 2 (12; A+cr) in accordance with Theorem 12. 11 1 for a — /3 — 1q. 

Set A := 12 . Note that A C M iVd is connected and compact. In the following we use the product 
measure Idili Fi on + where p, := Aj + Uj is defined on 12 and the index i gives reference to the 
corresponding coordinate. For functions /', g G (7 1 (A), i G I and x J G 12 for j G /, j i, define 




,x N ) := G{f(x 1 , 


2-1 # 2+1 
1 ^ J 1 ^ ? 


,x N ),g(x 1 ,...,x 


2 — 1 


,x i+1 ,...,x N )). 


Define the symmetric bilinear form (£,29) by 

N r 

Ff,g)-=y] L £‘(f,s)U d n for/.jeP-c^A). ( 3 . 1 ) 

,.1 Ja ‘ & 

Using the definition of the form Q yields 

, /• N N 

£(/,s0 = r / (V*/, Viflf) + 5 l A «.r(Vr,i/,V r> tf)) n dpj , (3.2) 

2 Ja i=1 j =1 

^- v - ' 

=:r (/,<?) 

where A* ,n := {x = (xi,..., xn ) G A| Xj G 12} and A* ,r := {x = (xi,..., xn) G A| x* G T}. I 11 
particular, A ?,n U A 2 ’ 1 = A for every i = 1,..., N. 

Condition 3.1. g G Z+A; Hi=i Fi)> Q > 0 n^i Fi~ a - e -- 
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3 The Dirichlet form and the associated Markov process 


Define /i by /j := p nf=i Tj and (£7 ^) by 


£(f,g) ■= 


r N N 

/ E + 5 WV ri i/,Vr,tf)) pjjdh 

i=l j=l 


(3.3) 


1 

2 


AT 


r(/,p) gY[dfij 

j =i 


= 2 / ■’(/.//) d T 

Note that the case 5 = 0 correpsonds to the setting of a system of particles which has a sticky but 
static boundary behavior. Then, the bilinear form (£,T>) can be written in the simpler form 

1 r N N 

£(f,g) = x / E 1 a*.o (Vj/, v 4 p) e n for /’ 9 e z>. 

Z 'A i=1 J= l 

By the fact that p is a Baire measure on A we get the following result: 

Proposition 3.2. Under Condition \3.1\ we have that C°°( A) is dense in L 2 ( A; //). 

Dehne A# := {x G A| X; G D for i G B, Xi G T for i G /\ 5 } and z/# := Hies n i & i\B a i- Tli en 


A = U Bc/ Ab an d A=E 


e 


B 


In this terms it holds 


where 


£{f,g)= E S s(f,g) for f,gev, 

<6^BCI 


£b(L g)-=\ I 5Z(Vi/> v ^) + 5 E Vr^) d/IB. 

“'As i£l\B 

Moreover, dehne for x G T JV ~I S I, 5^0, and p G A 1 (A; n^i Ah) 

:= {*/ e n |B| 


p 1 TT Xi < oo for some e > 0 }. 

J {z&Q. I s l| |z-j/|<e} i&B 

The dependence of x is given in the sense that the variables of p given by the index set I\B 
are fixed by the components of x. Since p is an element of L 1 (A B \ y B ), R^(B,x) is only defined 
Uiei\B<U almost everywhere. Similarly, for y G let R^(B,y) be given by 


R^(B,y) := {x G T^ f p 1 TT cr* < oo for some e > 0}. 

J{zer»-\B\\ \z—x\<e} 

In this case, the variables of p given by the index set B are fixed by the components of y and 
R l g (B,y) is only defined A* almost everywhere. Note that in both cases B determines the 
components which are not at the boundary. 

The following condition is a generalized version of the usual Hamza condition (see e.g. |MR92l 
Chapter II, (2.4)]): 
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Condition 3.3 (Hamza condition). It holds 

(HI) g = 0 riie-B onfl^\R^( K B,x)forYl i£l ^ B a i -a.e. x E r iV_ l B l for every 0 ^ B C / 

and if 5 = 1 additionally 

(H2) g = 0 El i eI \B a i- a - e - on r N ~^\R r Q (B,y) for Y[, &B \i-a.e. y E ff |Bi for every BEI. 

(For B = I the condition (HI) and for B = 0 the condition (H2) reduce to the ordinary Hamza 
condition.) 

Remark 3.4. (i) Condition 13.31 is a natrual generalizaion of the ordinary Hamza condition, since 

in the present setting of sticky particles we are also interested in dynamics whenever one 
(or several) particles are located at the boundary. The set B determines the components 
inside Q and its complement I\B the components on T. Tims, (HI) ensures that the Hamza 
condition for the components inside hi is fulfilled, wherever the remaining components stick 
on T. Since we are also interested in dynamics on T if <5 = 1, (H2) is the corresponding 
condition in this case. 

(ii) For ff = (0, oo) (HI) of Condition 13.31 coincides with (FGV14i Condition 2.7] (disregarding 
that (0, oo) is unbounded), since in this case the surface measure on T reduces to the case of 
the point measure in 0. 

Remark 3.5. If g is e.g. continuous on A and positive ]Q^ =1 /i*-a.e., then g is outside the set 
{£> = 0} locally bounded away from zero and hence, R^(B,x) = {y E g(z(B,x, y )) > 0} and 
R F e (B,y) = {x E g(z {B ,x, y )) > 0}, where 

f y' ll{i) , if iEB 

LB,x,y) | x i-7fl(0 ) if i E I\B 
with : I —> {1,..., \B\}, i (->■ |{1 < j < i\ j E B}\. Hence, 

\R*{B,x) = {y E L! |B I| g(z {B , x ,y)) = 0} 


and 

r N ~W\lf e (B,x) = {y E r iV -l B l| g(z {B , x , y) ) = 0} 
for every x E r JV ~l' B l, y E and Condition 13.31 is fulfilled. 

Lemma 3.6. Suppose that Condition 13.11 and Condition 13.31 are satisfied. Then the bilinear form 
(£,T>) is closable on L 2 (A; /i). 

Proof. Let ( fk)keN be an £-Cauchy sequence in T> such that fk —?• 0 in L 2 (A; y) as k —> oo. In 
particular, (/ fc ) fce n is ^n-Cauchy and converges to 0 in L 2 ( A B ] ys) f° r every 0 ^ B C I. Thus, by 
definition of £ B we have that (djfk)ke n is Cauchy in L 2 (A B Lb) for every j = d{i — 1) + l, where 
iEB and l E {1,..., d} and hence, djfk —>• hj E L 2 (Ab ; y b ) as k —> oo. In other words, 


/rJv-lsi j^|s| 


(d-jfk ~h 3 ) 2 g\\ \ JJ (Ti 

i&B i£l\B 


—y 0 as k —» 


oo. 


(3.4) 


Therefore, it exists a subsequence (dj/fcJzeN such that djfk t —* hj as l —> oo in L 2 (V( B I; gY\i<zB 
n i6J \ s ^-a.e. and similarly, fk t —>• 0 as l —> oo in L 2 (ffl B l; ELe/ye <T-a.e.. This implies 
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that hj = 0 on £7^ oWi&B -Va.e. E[ iei\B (J ^ a - e - by (HI) of Condition 13.31 (see |MR92l Chapter II, 
Section 2a)]) and hence, hj = 0 pis- a.e. on A b- In the case 5 = 1, we obtain a similar statement 
for the components of Vr ,ifk, i G I\B, k E N, by considering the term in analogy to (13.4(1 and 
integrating first with respect to nie/\B and afterwards with respect to n*e_B 'V By Fatou’s 
lemma holds 

£ ( fk, fk ) < hm inf £ (f k - f k ,, fi, ~ fk d ) 0 as k -E oo. 

I—too 

□ 


We denote the closure of (£,D) on L 2 (A;/i) by (£,D(£)). 

Proposition 3.7. Suppose that Condition \3.1\ and Condition \3.3\ are satisfied. Then (£,D(£)) is 
a symmetric, regular Dirichlet form. 

Proof. By Proposition 13.21 and Lemma, [3.61 (£. T>) is symmetric, densely dehned and closable with 
closure (£,D(£)) which is also symmetric. Moreover, by |MR92l Chapter I, Prop. 4.10] and 
the representation (13.31) . (£,D(£)) possesses the Markov property. Finally, C' 00 (A) C C 1 (A) C 

D{£) n C( A) implies that D(£) D C( A) is dense in D(£) with respect to the £f -norm as well as in 
C(E) with respect to the sup-norm. Hence, (£,D(£)) is regular. □ 

Proposition 3.8. Suppose that Condition \3.1\ and Condition \3.3\ are satisfied. Then the symmet¬ 
ric, regular Dirichlet form (£,D(£)) is strongly local and recurrent. 

Proof. Using (FOTlll Theo. 3.1.1] and (FOT111 Exercise 3.1.1] it is sufficient to show the strong 
local property for elements in D. Therefore, let f,g E D such that g is constant on some open 
neighborhood U of supp(/) (in the trace topology of A). Then it follows immediately by (13.31) that 
£{f,g) = 0. Hence, (£, D{£)) is strongly local. Clearly, 1a E V El D{£) and £(1a, 1a) = 0. Thus, 
(£,D(£)) is also recurrent. □ 

We summarize the preceding results in the following theorem: 

Theorem 3.9. Assume that Condition Id. II and Condition 1 3.3\ are fulfilled. Then the symmetric 
and positive definite bilinear form (£,V) is densely defined and closable on L 2 (A;/r). Its closure 
(£,D(£)) is a recurrent, strongly local, regular, symmetric Dirichlet form on L 2 (A;/r). 

By the theory of Dirichlet forms, we obtain immediately the existence of an associated diffusion 
process. For details see e.g. |MR92I Chap. V, Theorem 1.11] or (FOTllI Theorem 7.2.2 and 
Exercise 4.5.1]. We remark that the definitions of capacities (and hence, of exceptional sets) used 
in the textbooks |FOTll] and |MR92| are introduced in different ways, but that the defintions 
coincide in our setting (see |MR92l Chap. III, Remark 2.9 and Exercise 2.10]). (T t ) t>0 denotes the 
sub-Markovian strongly continuous contraction semigroup on L 2 (A; p) corresponding to (£,D(£)). 

Theorem 3.10. Suppose that Condition \3.1\ and Condition \3.3\ are satisfied. Then there exists 
a conservative diffusion process (i.e. a strong Markov process with continuous sample paths and 
infinite life time) 

M := (n,E, (Ji)^o, (X t )t> 0 ) (®t)t> 0) (Px)o;Ga) 

with state space A which is properly associated with (£,D(£)), i.e., for all (p-versions of) f E 
B b (A) C A 2 (A; p) and all t > 0 the function 

A9m ptf(x) := E x (f{X t )) := / f(X t )dP x E R 

J A 
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3 The Dirichlet form and the associated Markov process 


is a quasi continuous version of T t f . M is up to p-equivalence unique. In particular, M is p- 
symmetric (p is stationary), i.e., 


p t f g dp = / / p t g dp for all f,gE Bb( A) and all t > 0, 


and has p as invariant measure (p is reversible), i.e., 


p t f dp = / f dp for all f G £^(A) and all t > 0. 


Remark 3.11. Note that M is canonical, i.e., ft = C(M + , A) and X t (u;) = oj E ft. For each 
t > 0 we denote by © t : ft — y 12 the shift operator defined by ©t(cv) — u j(- + t) for oj £ ft such 
that X s o © f = X s+ i for all s > 0. We take into account to extend the setting to C'(M + , WL Nd ) by 
neglecting paths leaving A. 


3.2 Densities with product structure 

We introduce a special case of the setting given in Section HOI which will be of particular importance 
later on. 

Condition 3.12. Assume that g is of the form 

N 

g(x) = 4>(x) Qi(x l ) for x = (x 1 ,..., x N ) £ A. (3.5) 

i —1 

Qi G L 1 (12; A + a), i — 1,..., N, is given as in \2.3\) for some a, G L 1 (f2; A), > 0 X-a.e. and 

/3i G L l (Y',a), Pi > 0 a-a.e. such that the respective Hamza conditions are fulfilled (see Section 
Moreover, is a n^Li(cbAj + PiCrp-a.e. positive, real valued, measurable function on A such 
that <fi G A 1 (A; + Picrf)). Furthermore, we assume that <f fulfills Condition \3.S\ . 

Remark 3.13. Note that Condition 13. 121 implies Condition 13.11 and Condition 13.31 

Under these conditions it is also possible to consider the form defined in (13.31) from a different 
point of view. Define the form (£*, D(£ 1 )), i = 1,..., N, as the closure of the bilinear form 

\ J(Vf, Vs) OidX + 6 - J{Vrf, Vr g) Pda for /, g G C 1 ^) (3.6) 

on L 2 {ft ; cpA + Pia) and set p t := cti A, + Pi(Ti. Then, it is possible to define (A, T>) and (A, V) as in 
(13.21) and (13.31) respectively with g replaced by <f>. This construction yields the same bilinear form 
(£,V) on L 2 (A; p ), where p = Q]jf =l (Xi + af). 

Roughly speaking, the first definition of (£,T>) in Section [All corresponds to a Girsanov transfor¬ 
mation of N independent sticky Brownian motions on Q with constant stickyness along T (each 
associated to the form (Q,D(Q))) such that the transformed process has a drift given by V In o. 
In the present section, the form (£ l , D (£*)), i = 1 describes a distorted sticky Brownian 

motion on D with drift V In a, inside D and the stickyness along T is given by ^ as well as a drift 
along T given by Vr In Pi. Then, the Girsanov transformation by <f> yields an additional drift V In p. 
Note that the resulting form and process (up to equivalence) are the same, since the pre-Dirichlet 
forms on C 1 (A) coincide. 
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4 Analysis of the Markov process 


Densities with product structure as presented in the present section have the advantage that we can 
handle the densities gi, i — 1,..., N, by considering the forms (£*, D(£ 1 )), i = 1,..., N, as given 


in (13.611 . For this type of Dirichlet form it is possible to use the (regularity) results of |GV14a| . 


In this way the assumptions imposed on g t , i = 1,... ,1V, are not very restrictive. Only for the 
interaction part 0 it is necessary to demand stronger requirements. 

4 Analysis of the Markov process 

4.1 Generators and boundary conditions 

By Friedrichs representation theorem we have the existence of a unique self-adjoint generator 
(L,D(L)) corresponding to (£,D(£)). 

Proposition 4.1. Suppose that Condition \3.1\ and Condition \3.3\ are satisfied. Then there exists 
a unique, self-adjoint, linear operator ( L,D(L )) on L 2 (A]g) such that 


D(L) C D(£) and £(f,g ) = (-Lf,g) L 2 ( A;/i) for all f G D(L ), g G D(£). 


In order to determine (L,D(L)) for a suitable class of functions we need the following additional 
condition on g and T: 


Condition 4.2. Assume that F is C 2 -smooth, g is given as in Section [3.31 Moreover, it holds 
0 G C 1 (A) such that V In 0 G L 2 ( A; /i). cq,/3j G (7(11), yjofi G // 1,2 (D) and if 5 — 1 \ffii G H 1,2 ( r) 


for i — 1,..., N. 

Remark 4.3. Note that if cq, i = 1,..., N, and 0 are a.e. positive and the additional conditions 


of Condition 14.21 arc fulfilled, Condition 13.11 and Condition 13.31 are implied in view of Remark 13.51 


Proposition 4.4. Suppose that Condition\4.3\ is satisfied and let f G C 2 (A). Then 


N 



(4.1) 


where L l,Q f, L l,r f and L 1 ^ f for i = 1,..., N are given by 



(4.2) 


(4.3) 


(4.4) 


(4.5) 


Proof. Let / G C 2 ( A) and g G X? = C 1 (A). By integration by parts and (13.51) follows 



N N 


i =1 3 =1 
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4 Analysis of the Markov process 


\Y J i (VJ, V i9 ) gf[(d\ j + daj) + 5 -Y j I (Vr ,if, Vr ,g) efftdA,- + da 3 ) 

2=1 J = i 2=1 J = 1 

N 

5 E L_, ( / < v '/- v < 9 ) « "<) II(^ + 

Z i=l An An j/j 

S N f f 

+ 2 x /-iv-i (/ ( Vr >*/’ Vr -^) ^ nc^ + dajS) 

i =1 An Ar 

- [ (Aj/ + (—, V*/))p gdAi + f ^ (ra, Vj/)(? £>dcq) TT(dAj + her. 
An £ Ar Pi 7 ^ 


^X 


I /.*_i (- /( a ^/ + Vr ’*^) ^ rfa *) n^- +da - 

i=i ,7n Ar Q 


Note that for x = (x 1 ,. • ■ ; x^) £ A with x* £ 12 holds ^ = Ah- 1 + and similarly, if 6 = 1 and 
x* £ T holds ^Aii! = Vr V /3i + due to (13.5p . Hence, 


N 


£(/, g) = / - E ( 1 A , -° ( ii,n / + L ft ) + »*«• (V’ r / + £j r /)) 9 <*/l = / -Lf g dp 


2=1 


and therefore, the assertion holds true, since T> is dense in D{8). 


□ 


Remark 4.5. In contrast to the case of reflection or absorption, in the present case it is not necessary 
to require a boundary condition in order to determine the generator, since the reference measure 
contains the involved surface measures. Nevertheless, it is possible to derive a suitable condition 
such that the boundary terms vanish. More precisely, if we assume for simplicity that 5 = 0 and 
that / £ (A 2 (A) is given such that 


A</ + (—,V i /) + ^(n,V/) = 0 on A*’ r , i = l,...,iV, 
Q Pi 


(4.6) 


it holds 


£(/, S) = E / (A,;/ + (—, V<f)gdp - A f (A/ + (—, Vf))gdp, 

2 Ja Q 2 J A g 

i.e., the generator is of the well-known form. The condition (14.61) is called Wentzell boundary 
condition. Note that the drift in normal direction increases if the factor increases. Hence, it 
is justifiable to say that the boundary is less sticky for the 2-th particle at a point x £ T if /3j(x) 
decreases. This property can also be discovered in a similar way by [GV14al Corollary 4.17], since 
as a consequence of this ergodicity theorem the particle spends less time on the boundary if f r (dido 
decreases (compare also to |FGV14l Corollary 5.7]). Moreover, if we rewrite (14.61) in the form 

AA if + A(—V,;/) + ai(n, V/) = 0 on A*’ r , 

g 

and set A — 0 we obtain the Neumann boundary condition. 
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4 Analysis of the Markov process 


Define for i — 1,..., N 
as well as 


Ai \— 1L \i,nE -|- 5 1 \i,rP 


1/, /V,;a 7 ; V,;<?A _ t a.,, , Vr.iA , r Vr,i<, 


^ (~f l n + S 

where E denotes the d x d identity matrix. Then, set 


A := 


Using this notation, we get for / G C 2 ( A) the representation 

Lf — Ur(/lV 2 /) + (6, V/). 
Note that AA t = A 2 = A and in particular, P l = P. 

4.2 Solution to the martingale problem and SDE 


+ h 


■) . 


(Ai 

0 

... 0 \ 


(h\ 

0 

A-2 

0 

... 0 

and b : = 

1° 

0 A N j 


\bn) 


(4.7) 


(4.8) 


Theorem 4.6. The diffusion process M from Theorem, \3.10\ is up to p-equivalence the unique dif¬ 
fusion process having fi as symmetrizing measure and solving the martingale problem for (L, D{L)), 
i.e., for all g G D{L) 


g(X t )-g(Xo)- [ (Lg)(X s )ds, t> 0, 


is an T t -martingale under P x for quasi all x G A. Here g denotes a quasi-continuous version of g 
(for the definition of quasi-continuity see e.g. JMB,92[ Chap. IV, Proposition 3.3]). 

Proof. See e.g. |AR951 Theorem 3.4 (i)]. □ 

By Proposition 14.41 L is explicitly known on the set C 2 (A). Using the representation given in (14.81) . 
we obtain the following corollary: 

Corollary 4.7. Assume that Condition 4-2 is fulfilled. Let g G C 2 (A) and let M be the diffusion 
process from Theorem, \3. 1(A Then 

s(X t ) - g(X 0 ) - j^Tr(A(X s )V 2 g(X s )) + (b(X 8 ),Vg(X 8 ))ds, t > 0, 


is an IF t -martingale under P x for quasi every x G A, where A and b are defined in (4-3). 

Lemma 4.8 (weak solutions and martingale problems). Fix the probability measure P = x G A, 
on C'(M + ,M Ard ) (see also Remark [3.lip . Let A, b be given on A by (14.71) . If 

/(X,) - /(X„) - f 1 tV(X(X»)V 2 /(X,)) + (MX,), V/(X,)) ds 

is an J r /-niartingale under P for every / G Cf°(M. d ), the equation 

dX t = A(X t )dB t + b(X t )dt 

has a weak solution with distribution P, where (B t ) t >o is an Ahl-dimensional standard Brownian 
motion. 
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4 Analysis of the Markov process 


Proof. See e.g. |Kal971 Theorem 18.7]. Note that A and b (defined as in (14.71) ) fulfill the required 
conditions, since they are progressive. Furthermore, the property AA t = A is used. □ 

Remark 4.9. The solution to the SDE given in Lemma 14.81 results from M by extending the 
underlying filtration (J r t)t >o if necessary (see proof of |Kal97[ Theorem 18.7] and the references 
therein). For convenience, we use for the process equipped with the enlarged filtration again the 
notation 

M = (O, T , (^) t > 0 , (X t ) t > 0 i (O t )t>0, (Pi)i£a) 

taking into account that the associated Dirichlet form is still given by (£, D(£)). Since (X f ) t > 0 is 
A-valued, we use the notation 

X t = (X t \...,Xf), t> 0, 
where XJ is hi- valued for i — 1 ,,N. 


Theorem 4.10. M is a solution to the SDE 

dX‘ t =i„(x;)(dB; + i(Ap(x;) + x,))dt) - i r (xj)|i(xp n(x;)dt 

+Sl r (X;)(dB[- i +(^^(Xi) + ^(X,))dt), i=l,...,JV (4.9) 

dflp = P(Xj) o dBl 
Xo =x, 

for quasi every starting point i 6 A, where {B t ) t > 0; B t = [B ],..., B^), is an Nd-dimensional 
standard Brownian motion. 


Remark 4.11. A Fukushima decomposition of M (see jFOTlll Chap. 5]) yields the same result as 
in Theorem 14.101 We would like to mention that the argument used here in order to get a solution 
to the SDE (14.9ft does not work in this way for reflecting (Neumann) boundary conditions, since in 
this case the reflection is not given by a drift term. However, a Fukushima decomposition is still 
valid (see e.g. |Trnf)3| h because in this case it is also possible to assign an additive functional to 
the surface measure a. The advantage in our situation is that we are able to express the boundary 
behavior in terms of the generator. 


4.3 Solutions by Girsanov transformations 

Condition 4.12. For every i — 1,..., N, there exists Pi> 2 with pi > | such that 

- -— G L p , 1 (Vt D {Qi > 0}; cqA) and additionally —^- G LP (T D {& > 0}; /5jCr) if 5 = 1 

OLi Pi 

or equivalently 

^ e Ll c (Tt n { Qi > 0}; /q). 

OLi pi 

Moreover, cap £i ({gi = 0}) = 0. 

Define := > 0}. Assume that Condition 14.21 and Condition 14.121 are fulfilled. According 

to |GV14al Theorem 5.9] there exists for every i = 1,..., N a diffusion process 

M‘ := (PJ),>„, (Xj) t >0, (e;) t >0, (Pi) ie n,) 
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4 Analysis of the Markov process 


with strong Feller transition semigroup (p\)t> o and transition function (p l t (x, x £ D*. The 

processe M', i — 1,..., N, is associated to the form D(S 1 )) on A 2 (fb; nf), where /q = cqA + /3,cr. 
In particular, (p l t )t >o is absolutely continuous with respect to /q, i.e., for every t > 0 and x G fh, 
there exists a non-negative, measurable function p\(x,y), y G fh, such that 


P\{x-,A) = / Pl{x,y)dHi{y ) for every A G 13(0*). 


Let M be given by 


m := (Xf_, (x,) t20 , (e t ) t >„, (®",py*. 


=(x 1 ,...,rr-^)EA 




where A := x^Oj as well as 

X t (u) := (Xj^i),..., Xf (wjv)) and 0 t (w) := (©^i), • • •, ©f M) 

for u = (a>i,..., cujv) G Set P x := (gi^P^ for x = (x 1 ,... ,x N ) G A. 

Denote by (pt)t >o the transition semigroup and by (pt(x, -))t>o, x G A, the transition function of 
M. Then, it holds for every A = Ai x • • • x A N G C 13(A) 


Pt(x, A) = 




l A (Xi(w)) cIP x (uj) 


N 


' x N cyi . ", 

x ;=i“ i =1 

N 


[ l A (Xj(wi)) dP ,H 


V „ V 

[ / l A (X;(a Ji ))dPi,(a J( ) = nP;(^,A 

i=l i = l 


by dehnition of P x . Since x (1^13(0*) generates £>(A), it holds 


. v w 

Pt(x,A) = / U$(a;*,y*) JJd/x»(y*) for 
d A i=\ i= i 


every A G 13(A). 


As a consequence, pt(x, •), t > 0, x G A, is absolutely continuous with respect to Idili A* an d 

Ptf{x X ,...,x N ) =p? for every / G 13 6 (A), (4.10) 

where 

Plfix 1 , ...,x N ):= plf(x\ ..., x* -1 , ■, aA +1 ,..., x Ar )(x*) 
and the order of the pj, i = 1,..., N, is arbitrary. 


Consider the symmetric bilinear form on A 2 (A; \\ i=1 Pi) given by 


N r 

(«£i £‘)(f,g):=J2 

i=l 


dyj, 


where 


N 


/, g GD((8)^ 1 A 1 ) := {/ G A 2 (A; //j) | for each i — 1,..., N and for /q — 


a.e. 


i=l 
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4 Analysis of the Markov process 


(x 1 ,...,x t 1 ,x t+1 ,...,x N ) G : /(x\...,x® V,x* +i ,...,x n ) G D(£ 1 )}. 


N 


J+l 


Due to | BH91 . Chapter V, Section 2.1] {®hL l £\ D(®f =l £ 1 )) is a Dirichlet form on L 2 ( A; n^li Ah)- 
Obviously, this Dirichlet form extends the pre-Dirichlet form (£,T>) defined in (13.11) . 

Lemma 4.13. C' 1 (A) is dense in w.r.t. ((gi^f®)*, i.e., ((gjA^P, is the closure 

of (£,£>) on L 2 (A:Jlf =1 Hi). 

Proof. First, note that C' 1 (A) C D(®f =l £ l ) by dehnition of D(®f =l £ l ). For simplicity, we only 
consider the case N = 2. The statement for arbitrary IV G N follows by the same arguments. By 
|BH911 Proposition 2.1.3b)] D(£ 1 )®D(£ 2 ) is dense in D(£ l ®£ 2 ). Thus, it is sufficient to show that 
C 1 (A) is dense in D^ 1 ) <g) D(£ 2 ). Then the assertion follows by a diagonal sequence argument. 
Let h G D (P 1 ) <g) D(£ 2 ) such that ^(x^x 2 ) = f(x 1 )g(x 2 ) for n^ =1 /q-a.e. (x^x 2 ) G A, where 
/ G D(£( 1 ) and g G D(£ 2 ). Since C 1 (D) is dense in D(£ l ) and D(£ 2 ), we can choose sequences 
(fk)ken and ( gk)keN in C 1 (D) such that f k —> f in D(£ l ) and gk —>• g in D(£ 2 ) as k —> oo. Deffiie 
hk{x 1 ,x 2 ) := fkix^g^x 2 ) for x^x 2 G D. Then it follows easily by the prodcut structure of the 
underlying measure that the sequence (/i fc ) fceN , h k G C 1 (A), converges in A 2 (A; nLi Ah) ^ and 
moreover, the sequence is £ l ® £ 2 -Cauchy. □ 


Denote by (T t ®) t>0 the L 2 (Qj;/Xj)-semigroup of (£*,£)(£*)), i — 1,... ,N. By |BH911 Chapter V, 
Proposition 2.1.3] the L 2 (A; A 4 *)-semigroup (T t )t >o associated to {®^ =l £ l , D{®^ =1 £ 1 )) is given 

by 

N 

Ttf = If ■ ■ • T}f for/eL 2 (A; lift), 

2=1 

where 

■—T i f('r 1 P - 1 . P +1 v n \(P\ 

± t J \ x ) • • • > x ± t J P i ■ ■ ■ ■> x i M ) • • • j x J 

for x = (x 1 ,..., x n ) G A. Since M ! is associated to the form (A®, D(£ 1 )) for i — 1,..., N, it follows 
by (14.101) and Lemma [4.131 the following: 

Proposition 4.14. The Dirichlet form associated to M is given by the closure o/(A,C' 1 (A)) on 

L 2 ( A;nf=i^). 


Additionally to Condition 14.21 and Condition 14.121 we assume the following: 


Condition 4.15. <f is strictly positive. 


Under these conditions on (j) it is possible to perform a Girsanov transformation of M. Consider 
the multiplicative functional (Z t )t> o, Z t = exp (M t - ^), given by 

M t f \7 In <fi(X. t )dBt, t> 0. 

Jo 

Note that V In 0(X t ) = ^(X 4 ) and B tl t > 0, are M. Nd valued and also that V In 0 is bounded due 
to Condition 14.151 


In view of Remark 13.111 (applied to (£®, A>(£®)), i = 1,..., N), it holds x^jO* = C'(M + , A) and 
<8>ili-F = B(C(R+, A)). Thus, (x^O\ (g^jC) i s a standard measurable space (see |IW89l 
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4 Analysis of the Markov process 


Chapter I, Definition 3.3]) and hence, by [IW89 1 Chapter IV, Section 4] there exists for every 
i 6 A a probability measure Pf such that ( I®)!, L J-i = P i,t, where 


p UA) := / ^(wJdP^w) for A e x F t 


Let 


:= ( (X,) t >o. (0i)i>o, (P^Wt 


Then, the transition function ( pf(x , -))t>o °f is absolutely continuous with respect to p for every 
x G A. Indeed, by the previous considerations the transition function (pt(x,-))t >o is absolutely 
continuous with respect to p t . Assume that A e H(A) is given such that p(A) = 0. Since <f> is 
bounded from above and from below away from zero in view of Condition 14.151 and the continuity 
of (f, it also holds that flili /h(A) = 0 and hence, pt(x, A) = 0 for every t > 0 and x G A, i.e., 



1a (X t ) dP x 


0 for every t > 0 and x G A. 


Therefore, we also have 

pt(x,A)= I l A (Xt)dPi 



1a(x 4 ) dPi it = 



Zt 1 a(X*) dP x 


0 


We summarize the results of this section in the following theorem: 
Theorem 4.16. is a solution to the SDE 


dX‘ t =1 n (X*)(diJj + i(^i(Xj) + ^(xo)df) - 

+5 l r (X))(dS[’ J + (^i(Xj) + ^(X*))di 
dB V t l = P(Xj) o dB\ 


l r (X*)^(Xj) n(X*)di 


X 0 ^7 


(4.11) 


for every starting point x G A, where (B t ) t > 0? Eh = (D),..., is an Nd-dimensional standard 

Brownian motion. Moreover, the Dirichlet form associated to is given by (£, D(£)) on L 2 (A; /i) 
and its transition function ( pf(x , -))t>o is absolutely continuous with respect to p for every x G A. 


Proof. Due to results in |GV14a] every M ! solves the respective cZ-dimensional SDE for every 
starting point in Dj, i = l,...,iV. Hence, the process M solves the SDE for N independent 
particles, i.e., it solves (14. lip for (f given by the indicator function on A. As a consequence 
M'A solves (14. lip by the Girsanov transformation theorem (see |IW89L Chapter IV, Section 4]). 
Moreover, by the same arguments as in [GV14bj the Dirichlet form of the transformed process M c ' 
is given by (£, D(£)). □ 


4.4 Connection to random time changes 

In the following, we present the connections to random time changes for the case 5 = 0 and in 
particular, how the Dirichlet form construction is related to it. 
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4 Analysis of the Markov process 


As already mentioned in |GV14a| the sticky boundary behavior is strongly related to random time 
changes. Denote by Y = (Y t )t>o reflecting Browian motion on 14. The associated SDE is given by 

dYt = dB t + i dLj , 

where (Lf) t > 0 , Lf = — n(Y s )dlf , denotes the boimdary local time of Y. In this case, we 
also refer to this kind of boundary behavior as instantaneous reflection, since the process does not 
spend time on the boundary T. The underlying Dirichlet form is given by 


1 

2 Jn 


(V/, \7g) dX for f,ge hT’ 2 (14) 


on L 2 (14; A) and (If) t >o is the additive functional in Revuz correpsondence with the surface measure 
a on T. Note that a is in this case singular with respect to the reference measure A in the sense that 
the support of cr has measure zero with respect to A and this singularity describes the instantaneous 
reflection. Denote by (T)) t > 0 the inverse of the additive functional A t := t + /3(Y t )lf , t > 0, where 
f3 G (7(14) is strictly positive. Using ( T t )t>o as a new time scale it is possible to perform a random 
time change, namely we define the process X = (X t )t>o by Xj := Y T t for t > 0. Then, similar to 
the case of the positive half-line in |EP14| we obtain with the definition if- V£ t , t > 0, 


X t - X 0 = Y Tt - Y 0 = B Tt 
where Lf = — f* n(X s )dlf, and moreover, 

(B Tt ,B Tt ) = T t = f ln(Y. s )d, 

Jo 


\Ll = Br, + \ Lf, 


cTt 


S = 


ln(Y s )dA s 


= [ l n (Y T J ds= [ 1 q(X. 
Jo Jo 


)ds 


Thus, by eventually enlarging the filtered probability space there exists a standard Brownian 
motion (B t )t> o such that 

A i 

X 4 — X 0 = / 1 n (X s )dB s + -Lf. 


Furthermore, it holds 


^(X s 


lr(X s )ds — / 
Jo 


0(Y r ( 

T ‘ 1 


lr(Y Ts )dA Ts 


f -Y_lr(Y»)^ s = f l r (Y ,)dlj = % = if, 


i.e., /3(X. t )dlf = lr(X t )df and in particular, Lf = — lr(X s )n(X s )ds. As a consequence 

(X t ) t>0 solves the SDE 


dX t — lfi(X t )dB t — 


2 j3(X t 


l r (X t )n(X t )dt. 
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5 Application to particle systems with singular interactions 


According to the theory presented in [CFllj (see also |GV14aj ) the Dirichlet form corresponding 
to the time changed process is given by the closure of 


\ f_(Vf,Vg)d\ 



(V/,V^) l n (dX + /3da) 


for /, g G C' 1 (G) on L 2 (G; A + /3a). 


If we choose e.g. a 6 G 1 (fl) positive, a drift transformation by |Vlna yields a solution to 


dX t 


In (X- t )dB t 


1 

2 


l n (X t )Vln a(X t )dt 


1 a(Xf 

2^PQ) 


l r (X t )n(X t )dt, 


which is associated to the closure of 


F /(V/, Vg)°id\= ]- /(V/, Vg) ln(adX + /3da) for /, g e C 1 ^) on L 2 (G; aA + /3a). 

1 Jn 1 Jn 

Actually, the order of time change and Girsanov transformation does not matter in this case. This 
construction yields a single particle diffusing in G with sticky boundary behavior. Consequently, 
the idea in order to construct an interacting particle system with sticky boundary is to consider 
N independent particles in G which are connected to the tensor product of the forms (S l , D(S 1 )), 
i = 1 ,..., N, as presented above. Afterwards, a drift is introduced by the density </> which finally 
leads to the form considered in the present paper. 

An evident idea would also be to construct an interacting particle system with instantaneous 
reflection and to realize afterwards a time change. However, this is not possible in a simple way. 
The canonical Dirichlet form is given by the closure of 


\ l(Vf.Vg)t>d\ N for f,g 6 C‘(A) on L\A; e X N ), 


(4.12) 


where X N denotes the Lebesgue measure on A. For this kind of Dirichlet form we have a well- 
known regularity theory at hand which enables us to construct solutions to the underlying SDE 
even for singular drifts for every starting point in a specified set of admissible initial values (see 
e.g. |FG08| . |BG14 j and |FT95j ). Usually, only starting points in the corners of A (two or more 
particles at the boundary of G) are not admissible, since the boundary is not sufficiently smooth 
at these points. Nevertheless, such kind of dynamics do not diffuse on the boundary of A and 
hence, a time changed process will also not have this property. Therefore, it is not possible to 
construct an interacting particle system with sticky reflection via time change in use of the closure 
of 04.12j) . since a particle which reaches T is expected to sojourn a positive amount of time on T 
and meanwhile, the remaining particles keep on moving undelayed. This implies a diffusion on the 
boundary of A. An appropiate approach for a process with boundary diffusion and instantaneous 
reflection is given in [Tom80] and |Car09] . 

In |GV14a| it is shown that the transition semigroup of M‘, i = 1,... ,.N, given above has the 
strong Feller property. It seems not clear that the transition semigroup of M is doubly Feller 
(i.e., it is a Feller process with strong Feller transition semigroup). In this case, it would even be 
possible to deduce due to the results of | Chu86 j and |CK08] that the process 1VU of Theorem 14.161 
has the doubly Feller property. 


5 Application to particle systems with singular interactions 

In |Gra88| the author investigates a martingale problem with Wentzell boundary conditions in a 
very general form. In particular, the relation to SDEs is developed and an existence result is shown. 
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5 Application to particle systems with singular interactions 


As an application the author constructs a system of interacting particles in a domain with sticky 
boundary. This particle system gives a model for particles diffusing in a chromatography tube. 
More precisely, the considered domain is given by 0 := {x G R d | x\ > 0} and the investigated 
SDE on 0 reads as follows: 

dX t = a(X t )dN t + b(X t )(dt - p(X t )dK t ) + 'y(X t )dK t + r(X t )dC Kt , 

Xo = x G 0, 

where (X f ) t > 0 is a continuous, 0-valued process, (C t )t> o is a d- dimensional standard Brownian 
motion, (N t ) t > 0 is a d-dimensional continuous martingale and (K t ) t > o is given such that K 0 = 0, 
K t is increasing, dK t = lg@(X t )dK t , and 


(TV*, N j ) t = 8ij(t — [ P (X s )dK s ). 

Jo 

Here, the main focus is placed on the very general form of the martingale problem and SDE as well 
as the assumptions on a and a = oo T , which is not necessarily strictly elliptic. In former results 
(see e.g. [IW891 Chapter IV, Section 7]), it is assumed amongst other things that on > c > 0. In 
|Gra88| it is shown that the martingale problem with the sojourn condition p(X t )dK t < Iqq (Xt)dt 
has a solution if and only if the above SDE has a weak solution. Sufficient conditions are r = 0, 
a and b are uniformly Lipschitz continuous and bounded, 7 = n is the inward normal vector and 
p is bounded, measurable and positive. Nevertheless, the smoothness conditions on b are rather 
strong. If we assume additionally that an > 0 (e.g. if a is given by the identity matrix), it holds 
that 

p(X t )dK t = 1 de(Xt)dt. 

In the case of the identity matrix, the underlying SDE is given by 


dX t = l e {X t )dB t + b(X t )le(X t )dt + —— n(X t )dt, 

d( x <) 

Xq = x G 0 , 


where (B t )t> 0 is a d-dimensional standrad Brownian motion. This setting corresponds to the one 
considered in [ GV14a| for S = 0. The corresponding system of interacting particles is given by 

dX\ = 1 e(Xl)dBi + £7(X t )l e (Xj)df + —n(X*)dt, i = 1,..., N, 

P ( X H 

X 0 = x G 0' V , 

where X t = (X),..., X^). According to |Gra88] an application for this system of SDEs is a model 
for molecules diffusing in a chromatography tube. The particles are pushed by a flow of gas and 
are absorbed and released by a liquid state deposited on the boundary of the tube. Hence, it is 
resonable to suppose a sticky boundary behavior. However, it is physically unreasonable that two 
molecules are located at the same position in 0 at the same time. I 11 order to avoid this kind of 
behavior it is necessary to consider a singular drift b 1 , i — 1,..., N, which causes a strong repulsion 
if two particles get close to each other. The construction of such kind of stochastic dynamics via 
Dirichlet forms has already been realized for absorbing and reflecting boundary conditions. 


21 













5 Application to particle systems with singular interactions 


In analogy to |FG08 ;, Section 5], a continuous pair potential (without hard core) is a continuous 
function ( : W l —> M U {00} such that £(— x) = £(x) G M for all x G M d \{0}. £ is said to be 
repulsive if there exists a continuous decreasing function 77 : (0, 00) — y [0, 00) with lirn^o rj(t) = 00 
and R > 0 such that 

£(x) > v(.\ x \) f° r M < -R- 

In particular, £(0) = 00. For N E N and a repulsive continuous pair potential £ we consider the 
the function 

<j)(x) := exp(— ^ ((x l — x j )) for x = (x 1 ,..., x N ) G A = 0, N . 

l<i,j<N 

Note that <p{x) = 0 if there exist i. j G {1,..., N} such that x l = xK 

Let T be C' 2 -smooth. We assume that £ is a repulsive, continuous pair potential such that £> G C' 1 (A) 
and moreover, we assume that 

N 


V In cj) G L 2 (A; jj) with /! = (j> + A°i)> 


i =1 


where 07 and A are continuous and a.e. positive such that yTq G iL 1,2 (h2) and \/A £ iL 1,2 (r) for 
i — 1 ,..., N. Then Condition 14.21 is fulfilled and Theorem 14.101 can be applied, i.e., there exists a 
solution to the SDE 


dX] =Mxi)(dBi + l(^p(x;) - 5>,c(x; - x{))dt) - i r (x;)^(x;) n(x i t) it 

1 dAi 


+i Mxjjpsp + (Znh(xj) - V v rii c(x; - xf))dt 


A 




r A—1 
A 

* = 1,... ,1V 


= P(X‘) o dB,< 
X 0 =X, 


for quasi every starting point x G A. 

Example 5.1. A possible example is given by the Lennard-Jones potential 


C(*)=4e ((^) 12 -(o) 6 ), 

I A | 

where e and c are positive constants. It holds 


c 

Ixl 


V, In £>(x) = - V i C(x i - x J ) 
i¥=i 
24e 

> 121 ,. 

x* — x-q 


With /(r) := ^r(2 (^) U - (;) 8 ) we get 


V,ln0(x) = Y j /(|x® — x-^^x* — x j ). 

3+i 
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Thus, the absolute value of the acting force obviously depends only on the distance of the respective 
particles. In this case, for <5 = 0 the corresponding system of SDEs is given by 

dxi + * +W/(|x;-x||)(x;-x?) dt) 

1 j¥=i 

- lr(Xj)|(Xj) n(Xj) dt, i — 1,..., N, 

X 0 =x, 

where (B t )t> o, B t = [B],..., B^), is an Ahi-dimensional standard Brownian motion. It is natural 
that we obtain in this case only a solution for quasi every starting point, since points in A which 
describe configurations where two or more particles are at the same position in 0 are naturally 
not admissible in view of the singularity of C in 0. An appropriate regularity results regarding 
the elliptic PDE associated to the form (£,D(£)) would allow to apply the results of |BGS13j . 
In this case, a process on A\{</> = 0} = {x = (x 1 ,..., x N ) G A| x l ad for every i j} can be 
constructed which is a solution to the above SDE for every starting point in A\{</> = 0}. 
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